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Abstract 



We consider a linear integro-diffcrential equation which arises to describe both aggregation- 
fragmentation processes and cell division. We prove the existence of a solution (X,U,4>) to the 
related eigenproblem. Such eigenelements are useful to study the long time asymptotic behaviour 
of solutions as well as the steady states when the equation is coupled with an ODE. Our study 
concerns a non-constant transport term that can vanish at x = 0, since it seems to be relevant to 
describe some biological processes like proteins aggregation. Non lower-bounded transport terms 
bring difficulties to find a priori estimates. All the work of this paper is to solve this problem using 
weighted- norms. 

Keywords Aggregation-fragmentation equations, eigenproblem, size repartition, polymerization pro- 
cess, cell division, long-time asymptotic. 
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1 Introduction 

Competition between growth and fragmentation is a common phenomenon for a structured population. 
It arises for instance in a context of cell division (see, among many others, [Tj H"l IS] l6j [T4"l [T9 l [2Tj [25] [33] ) . 
polymerization (see CE]), telecommunication (see [2]) or neurosciences (see [30]). It is also a 
mechanism which rules the proliferation of prion's proteins (see |10[ [2"U1 I23j). These proteins are 
responsible of spongiform encephalopaties and appear in the form of aggregates in infected cells. Such 
polymers grow attaching non infectious monomers and converting them into infectious ones. On the 
other hand they increase their number by splitting. 
To describe such phenomena, we write the following integro-differential equation, 



0_ 

dt 




x^0 



< 




(1) 



[ u(0,t) = 0. 
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The function u(x, t) represents the quantity of individuals (cells or polymers) of structured variable 
(size, protein content...) x at time t. These individuals grow (i.e., polymers aggregate monomers, or 
cells increase by nutrient uptake for instance) with the rate t(x). Equation ([I]) also takes into account 
the fragmentation of a polymer (or the division of a cell) of size y into two smaller polymers of size 
x and y — x. This fragmentation occurs with a rate f3(y) and produce an aggregate of size x with the 
rate n(x,y). Equation ([1]) is a particular case of the more general one 

d d f°° 

—u(x,t) + —(T(x)u(x,t)) + [[3(x) + (i(x)]u(x,t)=nj j3(y)n(x,y) u(y,t) dy, x ^ x , (2) 

with the bound condition u(xo,t) = (see [TOl I23|). Here, polymers are broken in an average 
of n > 1 smaller ones by the fragmentation process, there is a death term fj,(x) > representing 
degradation, and a minimal size of polymers xq which can be positive. This more general model is 
biologically and mathematically relevant in the case of prion proliferation and is used in [91 \TU[ [20l [23] 
with a coupling to an ODE. Our results remain true for this generalization. 

A fundamental tool to study the asymptotic behaviour of the population when t — > oo is the existence 
of eigenelements (A, U, 4>) solution of the equation 



^-{t{x)U(x)) + (J3(x) + X)U(x) = 2 j°° P(y)K(x, y)U(y)dy, 
tU(x = 0) = 0, U(x)>0, f™U(x)dx = l, 



x ^ 0, 



-r(x)|-((Xx)) + (/?(x) + \)<f>(x) = 2/3(x) 



(3) 



(x) > 0, J °° <j){x)U{x)dx = 1. 

For the first equation (equation on IX) we are looking for D' solutions defined as follows : U € L 1 (M + ) 
is a V solution if Mip 6 C^°(M + ), 

i p oo poo . poo . 

T(x)U(x)d x ip(x) dx+ A / U(x)ip(x)dx= / (3(x)U(x)y2 / tp(y)n(y,x) dy — ip(x)\ dx. (4) 

Concerning the dual equation, we are looking for a solution (f> € W^^°(0,oo) such that the equality 
holds in 1^(0,00), i.e. almost everywhere. 

When such elements exist, the asymptotic growth rate for a solution to ([T|) is given by the first 
eigenvalue A and the asymptotic shape is given by the corresponding eigenfunction U. More precisely, 
it is proved for a constant fragmentation rate (3 that u(x,t)e~ xt converges exponentially fast to pU{x) 
where p = J uo(y)dy (see [221 For more general fragmentation rates, one can use the dual 

eigenfunction <j) and the so-called "General Relative Entropy" method introduced in |28[ 131] , It 
provides similar results but without the exponential convergence, namely that 

\u(y,t)e- xt - {u Q ,(t>)U(y)\(t>{y)dy 

1 1 t— >oo 

where {uQ,(j>) = f uo(y)(j)(y)dy (see [25| [25]). 

The eigenvalue problem can also be used in nonlinear cases, such as prion proliferation equations, 
where there is a quadratic coupling of Equation ([T]) or ([2J with a differential equation. In [51I1U [ [Td ] [5B] 
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for instance, the stability of steady states is investigated. The use of entropy methods in the case of 
nonlinear problems remains however a challenging and widely open field (see [34] for a recent review). 

Existence and uniqueness of eigenelements has already been proved for general fragmentation kernels 
K,(x,y) and fragmentation rates /3(x), but with very particular polymerization rates t(x), namely 
constant (r = 1 in [31]), homogeneous (r(x) = x^ in [26]) or with a compact support (Suppr = [0, xm] 
in [13]). 

The aim of this article is to consider more general r as |10tl37| suggest. Indeed, there is no biological 
justification to consider specific shapes of r in the case when x represents a size (mass or volume) or 
some structuring variable and not the age of a cell (even in this last case it is not so clear that ^tr = 1 , 
since biological clocks may exhibit time distorsions). For instance, for the prion proteins, the fact that 
the small aggregates are little infectious (see [2H [37] ) leads us to include the case of rates vanishing 
at x = 0. 

Considering fully general growth rates is thus indispensable to take into account biological or physical 
phenomena in their full diversity. The proof of [31] can be adapted for non constant rates but still 
positive and bounded (0 < m < t(x) < M). The paper [26] gives results for r(0) = 0, but for a very 
restricted class of shape for r. The paper [TJ] gives results for r with general shape in the case where 
there is also an age variable (integration in age then allows to recover Problem ([I])), but requires a 
compact support and regular parameters. Here we consider polymerization rates that can vanish at 
x = 0, with general shape and few regularity for the all parameters (r, (3 and k). 

From a mathematical viewpoint, relaxing as far as possible the assumptions on the rates t,k,/3, as 
we have done in this article, also leads to a better understanding of the intrinsic mechanisms driving 
the competition between growth and fragmentation. 

Theorem 1 (Existence and Uniqueness) Under assumptions (j5j)- (|13|) . there exists a unique so- 
lution (\,U,(f)) (in the sense we have defined before) to the eigenproblem ([3|) and we have 

A > 0, 

x a rU e L p (R + ), Va > -7, Vp G [1, 00], 
x a rU G W 1 ' 1 ^" 1 "), Va > 

3k > s.t. G L°°(1R + ), 

1 + x K 

r^GL^ c (lR + ). 

The end of this paper is devoted to define precisely the assumptions and prove this theorem. It is 
organized as follows : in Section [2] we describe the assumptions and give some examples of interesting 
parameters. In Section [3] we prove Theorem Q] using a priori bounds on weighted norms and then we 
give some consequences and perspectives in Section [H The proof of technical lemmas and theorem 
can be found in the Appendix. 
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2 Coefficients 
2.1 Assumptions 

For all y > 0, k(.,u) is a nonnegative measure with a support included in [0, y]. We define k on (M+) 2 
as follows : k(x, y) = for x > y. We assume that for all continuous function tp, the application 
fip-U 1 -^ f ^{x)K{x,y) dx is Lebesgue measurable. 

The natural assumptions on k (see [20] for the motivations) are that polymers can split only in two 
pieces which is taken into account by 

J k(x, y)dx = 1. (5) 
So n(y, .) is a probability measure and G L^ C (IR + ). The conservation of mass imposes 

xn(x,y)dx = ~, (6) 

a property that is automatically satisfied for a symetric fragmentation (i.e. K,(x,y) = n{y — x,y)) 
thanks to ([5]). For the more general model ([2]), assumption © becomes f xn(x,y)dx = ^ to preserv 
the mass conservation. 

We also assume that the second moment of k is less than the first one 

x 2 

K(x,y)dx < c < 1/2 (7) 

yZ 

(it becomes c < 1/n for model ©)• We refer to the Examples for an explanation of the physical 
meaning. 

For the polymerization and fragmentation rates r and /3, we introduce the set 

V := {/ > : 3fi, v > 0, limsupx _/i /(x) < oo and liminf x v f(x) > 0} 

and the space 

L\ :={/, 3a > 0, /eL^a)}. 

We consider 

/3 G L / 1 oc (]R + *) n P, 3a > s.l r G Lf oc (R+ x a °dx) n P (8) 

satisfying 

V.ff compact of (0, oo), Bm^ > s.t. r(x) > m# for a.e. x £ K (9) 
(if r is continuous, this assumption Q is nothing but saying that for all x > 0, r(x) > 0) and 

36 > 0, Suppp = [b,oo). (10) 

Assumption (|10p is necessary to prove uniqueness and existence for the adjoint problem. 

To avoid shattering (zero-size polymers formation, see [31 [23]), we assume 

3C>0,7>0 s.t. I K(z,y)dz < min(l,c(-) 7 ) and -^-eij (11) 
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which links implicitely r to re, and also 

- r e Ll (12) 

On the other hand, to avoid forming infinitely long polymers (gelation phenomenon, see I18|). we 
assume 

x(3(x) 

x—^+oo 



lim — -— = +00. (13) 
- t(x) 



Remark 1 In case when ([lip is satisfied for 7 > 0, then (|7|) is automatically fulfilled (see Lemma{3\ 
in the Appendix). 



2.2 Examples 

First we give some examples of coefficients which satisfy or not our previous assumptions. 
For the fragmentation kernel, we first check the assumptions ([5]) and ([6]). They are satisfied for 
autosimilar measures, namely k(x,u) = | K o(|)> with kq a probability measure on [0, 1], symmetric in 
1/2. Now we exhibit some Kq. 

General mitosis : a cell of size x divides in a cell of size rx and one of size (1 — r)x (see [27j ) 

i$> = \(5 r + 8 X - r ) ^ rG [0,1/2]. (14) 

Assumption fill jl is satisfied for any 7 > in the cases when r 6 (0, 1/2]. So (J7J) is also fulfilled thanks 
to Remark[TJ The particular value r = 1/2 leads to equal mitosis («(x, y) = 5 x= e). 
The case r = corresponds to the renewal equation (k(x, y) = ^(S x= o + S x=y )). In this case, we 
cannot strictly speak of mitosis because the size of the daughters are and x. It appears when x is 
the age of a cell and not the size. This particular case is precisely the one that we want to avoid with 
assumption (J7]) ; it can also be studied seperately with different tools (see [M] for instance). For such 
a fragmentation kernel, assumption (llip is satified only for 7 = 0, and the moments f z Ko(z)dz are 
equal to 1/2 for all k > 0, so (JT)) does not hold true. However, if we consider a convex combination of 
Kg with another kernel such as Kg with r £ (0, 1/2], then (]llj) remains false for any 7 > but ([7]) is 
fulfilled. Indeed we have for p € (0, 1) 

J z 2 (p4{z) + (1 - pK(z)) dz = P - + ^(r 2 + (1 - rf) = \{l - 2r(l - r)(l - p)) < \. 

Homogeneous fragmentation : 

K g'(^) = ^±l( z « + (1_^)«) for a>-l. (15) 

It gives another class of fragmentation kernels, namely in L 1 (unlike the mitosis case). The parameter 
7 = 1 + a > suits for (fTT|) and so ((7|) is fulfilled. It shows that our assumptions allow fragmentation 
at the ends of the polymers (called depolymerization, see [23], when a is close to —1) once it is not 
the extreme case of renewal equation. 

Uniform repartition (K,(x,y) = -io<x<y) corresponds to a = and is also included. 
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This last case of uniform repartition is useful because it provides us with explicit formulas for the 
eigenelements. For instance, we can consider the two following examples. 



First example : t(x) = tq, ft(x) = ftox. 

In this case, widely used by [20], the eigenelements exist and we have 



U{x) 



A -- 

o v 2 

<P{X): 



0^0, 



with X 



1 



(1 + X). 



Second example : t{x) = tqx. 

For such ft for which there exists eigenelements, we have 



-x. 



A 



and 



SyU(vY 



For instance when ft{x) = ftox n with n £ N*, then the eigenelements exist and we can compute U and 
<j) and we have the formulas in Table [TJ In this table we can notice that U(0) > but the boundary 
condition tU(0) = is fulfilled. 



n = 1 


A = r 


U{x) = f^e -°o x 


0(x) = 


n = 2 


A = r 


U(x) = J^e 2 -o x 

V ' y 7TT0 




n 


A = r 







Table 1: The example r(x) = tqx, ft{x) = ft^x 11 and uniform repartition k(x,u) = ^o<x<y 
gives the eigenelements solution to ([3]). 



The table 



Now we turn to non-existence cases. Let us consider constant fragmentation ft(x) = fto with an 
affine polymerization r(x) = tq + T\x, and any fragmentation kernel k which satisfies to assumptions 
d5J)- (JHJ) - We notice that (fT3|) is not satisfied and look at two instructive cases. 

First case : tq = 0. 

In this case assumption f)12[> does not hold true. Assume that there exists U € L 1 (IR + ) solution of 
([3]) with the estimates of Theorem [TJ Integrating the equation on U we obtain that A = fto, but 
multiplying the equation by x before integration we have that A = t\. We conclude that eigenelements 
cannot exist if t\ ^ fto- 

Moreover, if we take n(x,y) = -io< x <y, then a formal computation shows that any solution to the 

_ _Ma 

first equation of ([3|) belongs to the plan Vect{x ,x r i }. So, even if fto = t\, there does not exist an 
eigenvector in L . 



6 



Second case : tq > 0. 

In this case (I12p holds true but the same integrations than before lead to 

/xU(x)dx = — — - — . 
po - n 

So there cannot exist any eigenvector U € L l (xdx) for n > /5q- 



3 Proof of the main theorem 

The proof of Theorem [1] is divided as follows. We begin with a result concerning the positivity of 
the a priori existing eigenvectors (Lemma H]). We then define, in Section 13.21 a regularized and 
truncated problem for which we know that eigenelements exist (see the Appendix [B] for a proof using 
the Krein-Rutman theorem), and we choose it such that the related eigenvalue is positive (Lemma [2]). 
In Section 13.31 we gi ye a series of estimates that allow us to pass to the limit in the truncated problem 
and so prove the existence for the original eigenproblem (J3]). The positivity of the eigenvalue A and 
the uniqueness of the eigenelements are proved in the two last subsections. 

3.1 A preliminary lemma 

Before proving Theorem [H we give a preliminary lemma, useful to prove uniqueness of the eigenfunc- 
tions. 

Lemma 1 (Positivity) Consider hi and (f> solutions to the eigenproblem ([3]) . 

We define m := infjx : (x,y) € Supp /3(y)n(x,y)\ . Then we have, under assumptions (JSJ) , (J6j) , ((9j) 

x,y 

and ([TUP 

SuppU = [m, oo) and tU{x) > \/x > m, 
4>(x) > Vx > 0. 

// additionaly - € Lq, then 0(0) > 0. 

Remark 2 In case Supp k = {(x,y)/x < y}, then m = and Lemma{l\and Theorem^ can be proved 
without the connexity condition ()10p on the support of (5. 

Proof. Let xo > 0, we define F : x i-> r(x)U(x)e x o . We have that 

F'(x) = 2e J -o r ( s) ds / (3( y ) K ( x ,y)U(y)dy > 0. (16) 

So, as soon as tIA{x) once becomes positive, it remains positive for larger x. 

We define a := inf{x : t(x)U(x) > 0}. We first prove that a < |. For this we integrate the equation 
on [0, a] to obtain 

pa poo 

/ / P{y)K{x,y)U{y)dydx = 0, 

JO J a 
poo pa 

/ f3(y)U(y) / K(x,y)dxdy = 0. 

J a JO 
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Thus for almost every y > max(a, b), J Q a k(x, y) dx = 0. As a consequence we have 

[V If y 

k{x, y) dx = / k(x, y) dx < — I xk(x, y) dx = — 
J a a J ^ a 

thanks to (JSJ and ([B]), and this is possible only if b > 2a. 

Assume by contradiction that m < a, integrating ([3]) multiplied by cp, we have for all <p £ such 
that Supp ip C [0, a] 

r r 

(p(x)(3(y)n(x, y)U (y) dydx = 0. (17) 



By definition of m and using the fact that m < a, there exists (p, q) S (m, a) x (6, co) such that (p, q) € 
Supp f3(y)n(x,y). But we can choose 93 positive such that tp(p)U(q) > and this is a contradiction 
with (|17p . So we have m > a. 
To conclude we notice that on [0, m], U satisfies 

d x {T{x)U{x)) + \U{x) = 0. 

So, thanks to the condition t(0)U(0) = and the assumption ([9]), we have U = on [0, m], so m = a 
and the first statement is proved. 

For 4>, we define G(x) := (fi(x)e Jx o T < s ) . We have that 

G'{x) = -2e J *° T ^ p(x) / n(y,x)(f>(y)dy < 0, (18) 



so, as soon as (j) vanishes, it remains null. Therefore eft is positive on an interval (0, x\) with x\ € 
U {+00}. Assuming that x\ < +00 and using that x\ > a = m because J 4>{x)U{x)dx = 1, we can 
find X > x\ such that 

G (x)dx = -2 I I e Jx o T < s ) 4>(y)P(x)K(y,x) dy dx < 0. 



x\ J x\ JO 



This contradicts that 4>(x) = for x > x±, and we have proved that <j)(x) > for x > 0. 
If i € Lq, we can take xo = in the definition of G and so <f>(0) > or (f> = 0. The fact that (f> is 
positive ends the proof of the lemma. 
□ 



3.2 Truncated problem 

The proof of the theorem is based on uniform estimates on the solution to a truncated equation. Let 
rj, 8, R positive numbers and define 



ty, (a;) 



f] < X < 7] 

t{x) x > rj. 



Then t v is lower bounded on [0, R] thanks to ([9]) and we denote by fi = n(r],R) := infm m t X] . The 
existence of eigenelements (X^U^cj)^) for the following truncated problem when 5R < fi is standard 
(see Theorem [2] in the Appendix). 
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r a 



R 



— (r„(x)^(x)) + tf(x) + A*)W*(x) = 2 y /3(y)K(x,y)U*(y) dy, < x < R, 
t v U*(x = Q) = 8, U*(x)>0, fU*(x)dx = l, 







(19) 



^(x)^(x) + (/3(x) + A;)0j(x)-2/3(x) / K(y,x)<^(y)dy = 5(^(0), < x < R, 



dx 



k ^(i?) = 0, <^(x) > 0, / ^(x)W^(x)dx = 1. 

The proof of the theorem Q] requires A^ > 0. To enforce it, we take 8R = ^ and we consider R large 
enough to satisfy the following lemma. 



Lemma 2 Under assumptions (|5|), ([8]) and ()13p . £/iere exists a Rq > suc/i i/iat /or aZZ i? > i?o, */ 
we choose S = i/ien we /iaue A^ > 0. 



Proof. Assume by contradiction that R > and A£ < with <5 = Then, integrating between 
and x > 0, we obtain 



> A I U(y) dy 
Jo 



y + 2l f R ' p(y)K(z,y)U(y)dydz 



Jz 
R , i-x 



= 5 - t(x)U(x) + / P(y)U(y) dy + 2 U K (z, y) dz) /3(y)U(y) dy 

»/0 J x J 

> 5 - t{x)U{x) + f P(y)U(y) dy. 



Consequently 

r(x)U(x) >5+ I'M 
Jo r{y) 

and, thanks to GronwalFs lemma, 

t(x)U(x) > 8J° ^) dy . 



i~{y)U{y)dy 



But assumption (|13j) ensures that for all n > 0, there is a A > such that 

^-T-r > -> Vx > yl 
r(x) x 

and thus we have 

t(x)U(x) > 5x n , Vx > A. 
Thanks to (|8|) we can choose n and ^4 such that x~ n r(x) < j for x > A and then we have 

r-R r-R r-R n o 

1=/ U(x)dx> U(x)dx>5 ——dx>—(R-A) 
Jo Ja Ja r{x) R 

what is a contradiction as soon as R > 2A; so Lemma holds for Rq = 2 A. 

□ 
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3.3 Limit as 5 -»■ for U* and A J 

Fix 77 and let <5 — > (then i? — >• 00 since <5i? = §). 

First estimate: A^ upper bound. Integrating equation (|19p between and i?, we find 

A$<<5+ [ P(x)U*(x)dx, 



then the idea is to prove a uniform estimate on J pU^. For this we begin with bounding the higher 
moments f x a f3U^ for a > max (2, ao + 1) := m. 
Let a > m, according to ([7]) we have 



/" x a fx 2 1 

/ — k(x, v) dx < / — ttk(x, y)dx < c < -. 
J y a ~ J y 2 y ' ~ 2 

Multiplying the equation on by x a and then integrating on [0, R], we obtain for all A > 77 
x Q ((l-2c)/3(x))^(2;)(i2; < a f x a ~ 1 r n (x)U^(x)dx 



= a x a ~ 1 T v (x)U^ l (x)dx + a / x a-1 r(a:)Z<J(x) di 
< ai"- 1 " 00 sup {x ao T(x)} +u} A , a [ x a f3(x)U*(x)dx, 

where wa,q is a positive number chosen to have crr(x) < ujA,aXf3(x), Vx > A. Thanks to ([7|) and (fT3|) . 
we can choose large enough to have 0JA a ,a < 1 — 2c. Thus we find 

Va > m, 3,4 a : V77, <5 > 0, / x a /3(x)^(x) < n ' ; — := B a . (20) 

7 1 - 2c - UJ Aa ,a 

The next step is to prove the same estimates for < a < m and for this we first give a bound on 
t v U^. We fix p € (0, 1/2) and define x v > as the unique point such that j Q Xv -^^dy = p. It exists 
because (3 is nonnegative and locally integrable, and is positive. Thanks to assumption (|12p . we 



know that x n — > xq where xq > satisfies Jq° ^rnxdy = p, so x v is bounded by < x < x^ < 



Then, integrating (fT9j) between and x < x„, we find 



x. 



r^(x)^(x) < 5 + 2 f3{y)U d n {y)K(z,y)dydz 



< 5 + 2 J P(y)U°(y)dy 

i-Xrf roo 

= 5 + 2 P(y)l4(y)dy + 2 P(y)U*(y) dy 

JO Jx v 

< 5 + 2 sup {t v U*} H M- dy + ^ T y m f3(y)rf(y) dy 

(o,x v ) Jo T n(y) x ™ Jo 



< 5 + 2p sup {T v U 5 n } + —B m . 



(0,x v ) x r? 



10 



Consequently, if we consider 6 < 1 for instance, we obtain 

sup r„(s)Mj(s) < 1 + 2i W* m ;= c (21) 

as€(0,x) 1 — Zp 

so TjjZ/^ is uniformly bounded in a neighborhood of zero. 

Now we can prove a bound B a for x a f5lA^ in the case < a < m. Thanks to the estimates (I20p and 



(|2ip we have 

[ x a f3(x)U*(x)dx = [ x a (3(x)U^(x)dx+ [ x a p(x)U*(x) dx 

J Jo Jx 



< x a sup {tJ4} [ ^ldy + x 



— a— m i „TOo/„\7/5/ 



x m /3(x)Z^(x) dx 



(0,5?) 

< Cpx a + 5 m x a " m := B a . (22) 
Combining ([20j) and ([22]) we obtain 

Vq > 0, 35 Q : V??, 5 > 0, y x a /3(x)Z^(x) fix < B a , (23) 

and finally we bound 

X s v = 5 + J I3U*<5 + B . (24) 
So the family {A^},5 belong to a compact interval and we can extract a converging subsequence 



Second estimate : W 1,1 bound for x a T v U^, a > 0. We use the estimate (|23|) . First we give a L°°bound 
for 7>,Wi by integrating (fT9j) between and a; 

r^x^x) < 5 + 2 / /3(y)W*(y) < 6 + 2B := D . (25) 
J o 

Then we bound x a T r] U^ l in L 1 for a > — 1. Assumption (|13p ensures that there exists X > such that 
r(x) < xf3(x), Vx > X, so we have for R> X 

f-X pR 

x a T v {x)U 5 Jx)dx < sup{r^J} / x a dx + x a+l (3{x)U 5 Jx) dx 
(o,x) Jo Jx 



(0,X) 

sup 

(0,X) 



x a+l 

< sup {t v U%} — — + B a+1 :=C a . 
' a + 1 



Finally 

Va > -1, 3C a : V77, <5 > 0, y x a T ri {x)U 5 r) {x) dx < C a (26) 
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and we also have that x a U^ is bounded in L 1 because r € V (see assumption ©). 

A consequence of 421 and 426]) is that x%£<J is bound in L°° for all a > 0. We already have ([25 

and for a > 0, we multiply (fTUj) by x a , integrate on [0, x] and obtain 

/■if, />i? 
x a T v {x)Ufa) < a / y^r^U^y) dy + 2 y a (3(y)U^y) dy < aC a + 2B a := D a , 
Jo Jo 



that give immediately 



Va > 0, 3D a : V??, 5 > 0, supx%(x)W*(s) < D a . (27) 

x>0 



To conclude we use the fact that neither the parameters nor are negative and we find by the chain 
rule, for a > 

J^(x%(x)W*(x))|dx < a I x a - 1 T v (x)U*(x)dx+ j x a \d x (Tr,(x)U*(x))\dx 



<aj x a ~ L T v (x)U°(x)dx + X d v I x a U°(x)dx + 3 I x a /3(x)U°(x) dx (28) 
and all the terms in the right hand side are uniformly bounded thanks to the previous estimates. 



Since we have proved that the family {x a T v U^}s is bounded in VF 1 ' 1 (M + ) for all a > 0, then, because 
t v is positive and belongs to V, we can extract from {U^}s a subsequence which converges in L 1 (R + ) 
when 5 — > 0. Passing to the limit in equation (|19|) we find that 





— {T v {x)U n {x)) + (P(x) + \ v )U v {x) = 2 / P(y)K(x, y)U n {y) dy, 
W„(0) = 0, U v (x)>0, JU V = 1, 



(29) 



with \ T) > 0. 



3.4 Limit as r) — > for U rj and A,, 

All the estimates (|20p - (j28p remain true for 5 = 0. So we still know that the family {x a T ri lA r} } ri belongs 
to a compact set of L , but not necessarily {U^}^ because in the limit r can vanish at zero. We need 
one more estimate to study the limit 77 — ?• 0. 

Third estimate: L°° bound for x a r v U v , a > —7. We already know that x a T r) U rj is bounded for a > 0. 
So, to prove the bound, it only remains to prove that x -7 ^^ is bounded in a neighborhood of zero. 
Let define f^-.xt-t sup( 0iX ) T 4^n- ^ we integrate ([29]) between and x' < x, we find 

t v (x')U v (x') < 2^ J P{y)U 11 {y)K{z,y)dydz < 2^ J (3{y)U rj {y)K{z,y) dy dz 



and so for all x 



f v (x) < 2 j J f3(y)U v (y)K(z,y)dydz. 
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We consider x v and x defined in the first estimate and, using (jlip and (|12p . we have for all x < x r 
f v (x) < 2 ( [ f3(y)U v (y)n(z,y)dydz 



dy 



o 



2 / (3(y)U v (y) / n(z,y)dzdy 



o 

/■oo /r\7 

< 2 / /%)Z^(y)min(l ) C' / - 
Jo v 



= 2 [ X P(y)U v (y)dy + 2C H /%K(y) (-) 7 dy + 2C /%)Z^(y) (-) 7 dy 
We set Vr,(x) = x _7 / r) (x) and we obtain 

Z^" #(y) 

(1 - 2p)V r? (x) <K + 2C ^V a (y) dy. 

Jx T v {y) 

2Cp 

Ke 1 ~ 2 p 

Hence, using Gronwall's lemma, we find that V v (x) < — — - — and consequently 

2Cp 

J( e l-2p „ 

x _ %(x)i4,(a;) < := C, Vx € [0,x]. (30) 

1 

This last estimate allows us to bound U v by ^— which is in L\ by the assumption (jlip . Thanks to 
the second estimate, we also have that J x a lA^ is bounded in L 1 and so, thanks to the Dunford-Pettis 
theorem (see [8] for instance), {^jr? belong to a L 1 -weak compact set. Thus we can extract a subse- 
quence which converges L 1 — weak toward U. But for all e > 0, {x"^},, is bounded in VF 1 ' 1 ([e, oo)) 
for all a > 1 thanks to ([28)) and so the convergence is strong on [e, oo). Then we write 

/re roo 
\u n -u\ = J \U V -U\ + J \u n -u\ 



oo 



< 2c r 4^+ \U V -U\. 

Jo t(x) 



The first term on the right hand side is small for e small because ^ G ij and then the second term 

77—5-0 



is small for r\ small because of the strong convergence. Finally U v — > U strongly in L 1 (R + ) and U 



solution of the eigenproblem ([3]) . 
3.5 Limit as 8, 77 — > for 0^ 

We prove uniform estimates on <^ which are enough to pass to the limit and prove the result. 
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Fourth estimate : uniform <fP -bound on [0,^4]. Let A > 0, our first goal is to prove the existence of 
a constant C$ (A) such that 

Vt/,*, sup 0j < C (A). 

(0,A) 

We divide the equation on (jA by t v and we integrate between x and x v with < x < x v , where x v , 

bounded by x and x, is defined in the first estimate. Considering S < ^ 1 ~ 2p ^ (fulfilled for R > ^j: 

JL) 
2RJ 

~~T~\ [ K {z,y)(j) S Jz)dz + X v -(j> S J 

T v\V) Jo A* 



since 5 = we find 



2(l-2p) 



and we obtain 



$(x) < <j> d (x v ) + 2 ^ K(z,y)<f>°(z)dz + x v -<l>°{0) 

Jx T r)\V) JO A* 

< ^J(^) + sup {4> 5 V } (2 [ " 44 r y) dz + x v - 

(o,x v ) v Jo T v(y) Jo v 

SUp <f>*(x) < - 6*(Xr,). 

xe(o,x) l-2p-8x/n 1 

rx Z^ + A^ 
J a 



Using the decay of (j) rj {x)e - Tt ? , there exists C(A) such that 

sup <f>fa) < C(A)<f> s n (x v ). 

x€(Q,A) 



Noticing that / 4> s n {x)U^{x)dx = 1, we conclude 

1> / </> s (x)U%(x)dx > <j> s (x v ) / e h U S Jx)dx, 
Jo Jo 

so, as x v —7- xo and U(x)dx > (thanks to Lemma [1] and because xq > b > a), we have 

sup <f> s < C (A). (31) 

(0,A) 



Fi/t/i estimate : uniform ip^-bound on [A, 00). Following an idea introduced in [52] we notice that the 
equation in (|19p satisfied by is a transport equation and therefore satisfies the maximum principle 
(see Lemma H] in the Appendix). Therefore it remains to build a supersolution cf> that is positive at 
x = R, to conclude (j>^(x) < (f>(x) on [0, R]. 

This we cannot do on [0, R], but on a subinterval [Aq, R] only. So we begin with an auxiliary function 
Tp(x) = x k + 6 with A; and 6 positive numbers to be determined. We have to check that on L4o,-R] 

-T(x)—lp(x) + (\ s v + l3(x))lp(x)>2p(x) J K(y,x)<p(y)dy + 6<f>*(0), 

i.e. 

-kT(x)x k - 1 + (X 5 v + /3(x))7p(x)> (26 + 2 J K(y,x)y k dy)f3(x) + 50 5 V (O). 
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For k > 2, we know that f dy < c < 1/2 so it is sufficient to prove that there exists Aq > 

such that we have 

- kr(x)x k - 1 + (Aj + (3(x))(x k + 9)> (29 + 2cx fc )/?(x) + 5C (1) (32) 
for all x > ^4oi where Co is defined in (|3T|) . For this, dividing (|32|) by x fc_1 r(x), we say that if we have 

( i- 2t) !g> t+ *»:; c ;'"» (33) 

then (|32p holds true. Thanks to assumptions (|8j) and ()13|) we know that there exists k > such that 
for any 9 > 0, there exists ^4o > for which (|33|) is true on LAo, +oo). 
Then we conclude by choosing the supersolution cj>(x) = C °^^ Tp(x) so that 

t(x) > <f> s v (x) on[0,A ], 

and on [Aq ,R], we have 

-r(x)^(x) + (A* + P{x)j${x) > 2/3(x) /* K (y, xj$(y) dy + <^(0), 

> 0, 



(34) 



which is a supersolution to the equation satisfied by <j>\ * Therefore (j)^ < 4> uniformly in 77 and 5 and 
we get 

3fc, 9, C s.t. V77, <5, <j} s n (x) < (Cx k + 9). (35) 

Equation (|19p and the fact that 0^ is uniformly bounded in L^ C (M + ) give immediately that is 
uniformly bounded in L^ C (R + , r(x)dx), so in L^ c (0,oo) thanks to ©. 

Then we can extract a subsequence of {(j)^} which converges C°(0, 00) toward (f>. Now we check that 
<j) satisfied the adjoint equation of ([3]). We consider the terms of (|19p one after another. 
First (A* + P(x))4> s v (x) converges to (A + P(x))<j>{x) in L% c . 

For d x (j>^, we have an L°° bound on each compact of (0, 00). So it converges L°° — *weak toward d x (p. 
It remains the last term which we write, for all x > 0, 

I K(y,x)((f> 5 Jy) - <j>{yj)dy < \\<t& - ^\\ LaO f 0jX ) — 0. 
Jo 

The fact that J <j)U = 1 comes from the convergence L°° — L 1 when written as 

cj> s 7] (x)U%(x)dx = J ^^(l + x k )U 5 n {x) dx — ► y ^^(l + a^)W(x)dx = ^ tfW. 



At this stage we have found (X,U,4>) G M + x L 1 (M+) x C(R+) solution of ©. The estimates 
announced in Theorem [1] also follow from those uniform estimates. It remains to prove that A > 
and the uniqueness. 
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3.6 Proof of A > 

We prove a little bit more, namely that 

A > - sup{r(x)£Y(x)}. 

2 x >o 

We integrate the first equation of (J3j) between and x and find 



(36) 



< A / U{y) dy 
Jo 



-t{x)U(x)- / (3(y)U(y)dy + 2 

Jo Jo 

roo roo 




(3{y)K(z,y)U(y) dy dz 

noo 



-t{x)U{x) + 2 / p{y)U(v) dy 
Jo 

-t(x)U(x) + 2A, 



Hence 2A > t(x)U(x) and (f3"B"j) is proved. 



3.7 Uniqueness 

We follow the idea of [26]. Let (Ai,Wi, </>i) and (X2M2, $2) two solutions to the eigenproblem ([3]). First 
we have 

Ai Jux{x)<t> 2 {x)dx = J(-d x {T(x)U 1 (x))-/3(x)U 1 (x)+2 J p(y)K(x,y)U x (y)dy^<t) 2 {x)dx 



r{x)d x (j) 2 (x) - P(x)<fa{x) + 2P(x) J K(y,x)<f) 2 (y)dy)Ui(x)dx 
A 2 / Ui(x)4> 2 {x) dx 



and then Ai = A2 = A because J U\(j) 2 > thanks to Lemma [TJ 

For the eigenvectors we use the General Relative Entropy method introduced in |28t I29j. For C > 0, 
we test the equation on U\ against sgn(^ — C)<p\, 

> I 

^(x)-C t{x)U 2 {x)Mx) we find 



d x {r{x)Ui{x)) + (A + P(x))Ui(x) - 2 / /%) K (x,y)Wi(y) dy] sgn(^|(x) - C^Mx) dx. 



Deriving 



J d x (r(x)Ui(x))8gsx(^(x) -C^4>i(x)dx = J d x (^ ^(x) -C T(x)U 2 {x)^x{x) s j dx 



+ 

and then 



J d x (r(x)U 2 (x))^(x)sgn(^(x) -Cj^i(x)dx 



U 2 



x) — C d x (r(x)U 2 (x)(f>i(x)) dx 
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14 r f x 

jy{x)-C \ (3(x)K(y,x)U 2 {x)4>i(y)dy 



dx 



+2 



f3(y)K(x,y)U 2 {y)dy ^-(x)sgn(^(x) - Cj(j>i(x)dx 
- + /3(x))^(x)sgn(^(x) - C)u % {x)<h{x) dx, 



j d x {T(x)U 1 {x))sgn(^\x) -C^4>i{x)dx 



U 2 



U 2 (y)(j)i(x) dxdy 



+2 



P(y)K(x,y)U 2 (y)dy ^(x)sgn^(sc) - c)0i(x)dx 



(A + /3(x))^(x)sgn(^(x) - c)u 2 {x)Mx) dx. 



U 2 



M 2 



So 



= 2 



P{y)n{x,y) Yj-(y)~C - ^-(x)-C U 2 (y)<f>i(x) dxdy 
l U 2 U 2 J 



+2y y p(y)K(x,y)U 2 (y)dy — (x)sgn(^(x) - Cj0i(x)dx 



fi(y)K(x,y)Ui(y)dy sgn(^(x) - C 



sgn(-(x) 



C s 



l-sgn( ^(x)-C)sgnf^(y)-C') = on the support of k(x, y) for all C thus ^(x) = fy(y) 



Hence 

on the support of k(x, y) and 



(37) 
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We can prove in the same way that 4>± = 4>2 even if we can have IA = on [0, m] with m > 0. Indeed 
in this case we know that /3 = on [0, m] and so 



<t>i{x) = <fc(0)e J ° T(s) Vx G [0,m], i G {1,2}. 



4 Conclusion, Perspectives 

We have proved the existence and uniqueness of eigenelements for the aggregation-fragmentation 
equation dU with assumptions on the parameters as large as possible, in order to render out the 
widest variety of biological or physical models. It gives access to the asymptotic behaviour of the 
solution by the use of the General Relative Entropy principle. 

A following work is to study the dependency of the eigenvalue A on parameters r and (3 (see |27j). 
For instance, our assumptions allow r to vanish at zero, what is a necessary condition to ensure that 
A tends to zero when the fragmentation tends to infinity. Such results give precious information on 
the qualitative behaviour of the solution. 

Another possible extension of the present work is to prove existence of eigenelements in the case of 
time-periodic parameters, using the Floquet's theory, and then compare the new Xp with the time- 
independent one A (see [H]). Such studies can help to choose a right strategy in order to optimize, 
for instance, the total mass f xu(t,x)dx in the case of prion proliferation (see [10]) or on the contrary 
minimize the total population J u(t,x)dx in the case of cancer therapy (see 112]). 

Finally, this eigenvalue problem could be used to recover some of the equation parameters like r and 
(3 from the knowledge of the asymptotic profile of the solution, as introduced in jl5[ [35] in the case of 
symmetric division (r = 1 and k = 5 x= v), by the use of inverse problems techniques. The method of 
[35] has to be adapted to our general case, in order to model prion proliferation for instance, or yet 
to recover the aggregation rate r ; this is another direction for future research. 
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Appendix 



A Assumption on n. 

Lemma 3 Assumptions dS]),© and (fTTl) with 7 > imply that 

-(l-ri)y 

k(x, y) alx > 0, 



inf lim 

y r?— >0 



'/// 



which means that polymers undergo a decrease in the size during fragmentation process. As a conse- 
quence, assumption flTJ) holds true. 



Proof. With the first assumption ([5]) we have 



mi 



n{x,y)dx = / K(x,y)dx + 



0--ri)v 



>w 



K(x,y)dx+ I K(x,y)dx. 
(i-n)v 



The two other assumptions (jBJ) and (jlip allow to control the mass of K at the ends : 



k(x, y) dx 



1 



k(x, y) dx 



u 



> l-Crf 



'0—n)y 
1 [v x 

1 - v J(i- V )y y 
1 1 



k(x, y) dx 



k(x, y) dx 



2(1-7/) r^O 2' 



which gives the first assertion of the lemma. 
Now we can prove (|7|) : 



y x 2 



y 



k(x, y) dx < 



>w 



-k(x, y) dx + 



< 



y 

r (i-v)v x 

w y 



-k(x, y) dx 



x 

w y 

0--v)y 



(1- V )y y 

k(x, y) dx +(1 — 77) 
k(x, y) dx 



+ 



0--v)v x 2 



k(x, y) dx 



(i-v)y 



y 



-k(x, y) dx 



< 



2~ V 



k(x, y) dx. 



We use the first part of the proof to conclude. Taking 77 = min^i, ^ 4g ^i/ 7 j for instance, we obtain 



and the lemma is proved for c — 



0--v)y 
w 

2 — 48' 1=1 



K(x,y)dx > -, 
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B Krein-Rutman 



We prove existence of solution for the truncated equation (I19p . In this part r] and 5 are fixed (with 
5R < (J,), so we will omit these indices for r, A, IA and (ft but we keep in mind that t(x) > p, > 0. 
We use the Krein-Rutman theorem which requires working in the space of continuous functions (see 
|31j for instance). First we define regularized parameters as follows : 



T £ = p £ *T, /3 £ = p £ *f3, and Vy > 0, k £ (., y) = p £ * «(., y), 
where p £ (x) = -/?(§ ) with p G C£°((0, oo)), positive and such that J °° p = 1. Then we have the theorem 

Theorem 2 Under assumptions ([5])- (|13p on £/te parameters and for all e > 0, i/iere is a unique 
solution A E €R andU £ ,4> £ G C 1 ([0, i?]) to i/te regularized eigenproblem 



( d_ 

dx 



(t £ {x)U £ (x)) + (&(x) + A e )W e (x) = 2 / p e (y)K e (x, y)U £ {y) dy, < x < R, 

Jo 

r £ U £ (x = 0) = 8 f R U £ (y) dy, U £ (x) > 0, j R U £ {x)dx = 1, 



-r e (x)— ^(x) + (/3 e (x) + A e )^(x)-2/3 £ (x) y *s e (|/,x)0 e (y)dy = r e (O)^ e (O), < x < R, 

MR) = o, 4>e{x) > o, f*<f> £ (x)u £ (x)dx = 1. 

(38) 

Proof. We follow the proof of |31j . We define linear operators on E := C°([0, R]) to apply the 
Krein-Rutman theorem. 

Direct equation. For v > we consider the following equation on E 

— (n(x))H ——n{x)-2 ——K £ (x,y)n(y)dy = ——, < x < R, 

dx t £ {x) J T £ (y) t £ (x) ^ 

n(x = 0) = 5f R ^j- ) dy, 

and we prove that the linear operator A : / i— > n (solution of ()39p ) satisfies to the assumptions of the 
Krein-Rutman theorem. 

First step: construction of A. Fix / G E and for m £ E, we define n = T(m) G E as the (explicit) 
solution to 

— (n(x))H 7 ——n(x) = 2 -—-K £ (x,y)m(y)dy^ — , < x < R, 

ox t £ (x) J T £ (y) t £ {x) 

n(x = 0) = 5j R ^dy, 

We prove that T is a strict contraction. Therefore it has a unique fixed point thanks to the Banach- 
Picard theorem. This fixed point is a solution to 
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In order to prove that T is a strict contraction, we consider mi and mi two functions in E, we 
compute for n = n\ — rii, m = m\ — 1712, 



— (n(x))^ — — n(x)=2 ——K £ (x,y)m(y)dy, < x < R, 

ox tAx) J T £ (y) 



therefore 



7) 







^ = 0) = Sl R ^dy, 



Te(x) 



\n(x)\ < 2 



T s(y) 



K e (x,y)\m(y)\dy, < x < R, 



n(x = 0)<5f R \^§dy. 



After integration, we obtain 



|n(x)|e 



.in 



< s 



m(y)\ 



dy+ [ ef° re e I ^^-K £ (x',y)\m(y)\dydx' 

T e{y) Jo Jo T £ {y) 



and thus 



\n(x)\ < 5[ R ^dy + 

Jo T e(y) 



e J*' r e / ——K e {x,y)\m{y)\dydx 

Jo Te{y) 



1 

< \\m \\e- 



< \\m\\ E - 



R 



5R+ I e~ j *'— I p £ (y)K e (x',y)dydx' 



5R + 



R 



Pe{y)K e {;y)dy 



— 11— w ^-x') , I 



< \\m\ 



fi L 



5R + v 1 ||t £ ||i,oo 



Pe(y)n e (-,y)dy 



Because 5R < /j by assumption, we can choose v large so that k < 1 and we obtain 

\\n\\E < fc||m||£. 

Thus T is a strict contraction and we have proved the existence of a solution to (I39p . 



Second step: A is continuous. This relies on a general argument which in fact shows that the linear 
mapping A is Lipschitz continuous. Indeed, arguing as above 



|n(x)|e Jo T e 



and thus 



JO Te(y) JO 



Jo r e I @ £ y k £ (x' ,y)\n(y) \ dydx' + / e*' T e ^-7-7: 

JO Te(y) JO r e (x') 



'■ R I ffv)l 

\n(x)\ <k\\n\\ E + I hL \-7rdx' < k\\n\\ E + -\\j\\e- 



This indeed proves that 



T e(x>) 

R 



\n\\E < 



ll/lb- 
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Third step: A is strongly positive. For / > 0, the operator T of the first step maps m > to n > 0. 
Therefore the fixed point n is nonnegative. In other words n = A{f) > 0. If additionally / does not 
vanish, then n does not vanish either. Therefore n(0) = 5 J Q R ^Ardy > and thus 

nix) > n(0 +e Jo / e Jo ^ dx > 0. 

Jo r £ {x') 

Fourth step: A is compact. For < 1, the third step proves that n is bounded in E and thus 

d v + Pz f My) f 

— n = n+ / ——K{x,y)n(y)dy^ 

Ox r £ J T £ {y) t £ 

is also bounded in E. Therefore by the Ascoli-Arzela theorem the family n is relatively compact in E. 



Adjoint equation. A function is a solution to the adjoint equation of (|38|) if and only if <j)(x) := 
<f>{R — x) satisfies 







R 



T £ (x)—<f){x) + (P £ (x) + Xe)<f>(x) - 2/3 E (x) I K £ (y, R - x)<j> e (y) dy = 5<j) £ {R), < x < R, 



o 



I Mo) = o, 

(40) 

where f e (x) = t £ (R — x) and (3 £ (x) = j3 £ {R — x). Then the same method than for the direct equation 
give the result, namely the existence of A and (j) solution to (|40l ). 



Finally we have proved existence of (Xu,U £ ) and (A^,, (f> £ ) solution to the direct and adjoint equations 
of (I38p . It remains to prove that \y = A^ but it is nothing but integrating the direct equation against 
the ajoint eigenvector, what gives 



Aw / U £ 4> £ = \<t, I U £ 

□ 



j U £ (j> £ = A J t 



To have existence of solution for (|19p . it remains to do e — > 0. For this we can prove uniform bounds 
in L°° for U £ and <f> e because we are on the fixed compact [0,JR]. Then we can extract subsequences 
which converge L°° * —weak toward U and (j), solutions to (|19p because r £ and j3 £ converge in L 1 toward 
r and j3. Concerning k £ , we have that for all ip € J (p(x)K £ (x,y)dx — > J (p(x)n(x,y)dx Vy, and 
it is sufficient to pass to the limit in the equations. 



C Maximum principle 

Lemma 4 If there exists A$ > such that <j)> (f> on [0, Aq\ and <f> a supersolution of ([3|) on [Aq,R\ 
with 4>{R) > (j)(R), then ~4> > 4> on [0, R]. 
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Proof. The proof is based on the same tools than to prove uniqueness (see above) or to establish 
GRE principles (see (28j [29] for instance). 

We know that <p > <p on. [0, Aq] and that <\> is a supersolution to the equation satisfied by 4> on 
L4o )-/?]) i.e. there exists a function / > 8(f)(0) such that 

-T(x)p4>{x) + (A + P{x)j${x) = 2P(x) [ n(y, x]4>(y) dy + f(x), Vx e [A , R]. 



dx 

So we have for all x € [Aq, R] 

-r(x)^-((f)(x) - ~$(x)) + (A + P{x)){<t>{x) - = 2f3(x) £ n(y, x)((f)(y) - 0(y)) dy - f(x). 

Then, multiplying by H^ > ^, we obtain (see [31] for a justification) 

-r(x)^(0 - ^)+(x) + (A + /3(x))(0 - 4>) + {x) < 2f3(x) £ K (y, x)(<f> - ~4>) + {y) dy - f(x)t^(x), 

and this inequality is satisfied on [0,R] since ((f) — </>)+ = on [0,^4o]- 
If we test against U we have, using the fact that 4>(R) = < <p(R), 

f 11 ^ - 4>) + (x)^-(r(x)U(x)) dx + [ R (X + P(x))((f> - ^)+(x)U(x) dx 
Jo dx J 

pR pR pR 

< 2 (<j>-4>)+(y) (3(x)^(y,x)U(x)dxdy- f(x)t^(x)U(x)dx. 



But if we test the equation © satisfied by U against ((f) — <p)+, we find 

rR p, r R 

10 



/ ((f) - (f)) + (x)— (r(x)U(x))dx + / (A + f3(x))((f) - <p) + (x)U(x) dx 
Jo Jo 



-R 

(<t>~ 4>)+{y) / P(x)n(y,x)U(x)dxdy, 



and finally, substracting, 



f R 

°<-J Q f(x)l^(x)U(x)dx, 



so 



R 

t^(x)U(x) dx < 

and this can hold only if l^^; = or 0(0) = 0. But we deal with the truncated problem with 



t(x) > rj > 0, so - € Lq and (f)(0) > thanks to the lemma [TJ Thus = and the lemma H] 

proved. □ 



is 
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